Non-thermal leptogenesis with almost degenerate superheavy neutrinos 
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We present a model with minimal assumptions for non-thermal leptogenesis with almost degenerate 
superheavy right-handed neutrinos in a supersymmetric set up. In this scenario the gauge singlet 
inflaton is directly coupled to the right-handed (s)neutrinos with a mass heavier than the inflaton 
mass. This helps avoiding potential problems which can naturally arise otherwise. The inflaton 
decays to the Standard Model leptons and Higgs via off-shell right-handed (s)neutrinos and reheats 
the Universe. The same channel is also responsible for generating the lepton asymmetry thus 
requiring no stage of preheating in order to excite superheavy (s)neutrinos. The suppressed decay 
rate of the inflaton naturally leads to a sufhciently low reheat temperature, which in addition, 
prevents any wash out of the yielded asymmetry. We will particularly elaborate on important 
differences from leptogenesis with on-shell (s)neutrinos. It is shown that for nearly degenerate 
neutrinos a successful leptogenesis can be accommodated for a variety of inflationary models with a 
rather wide ranging inflationary scale. 



I. INTRODUCTION 

The consistency of the abundance of the light ele- 
ments synthesized during the Big Bang nucleosynthesis 
(BBN) requires that the baryon asymmetry of the Uni- 
verse (BAU), parameterized as tjb = {ub — n-o)/ s with s 
being the entropy density and ub is the number density 
of the baryons to be in the range (0.3 — 0.9) x 10~^° [Q . 
The asymmetry can be produced from the baryon sym- 
metric Universe provided three conditions are simultane- 
ously met; B and/or L- violation, C- and CP- violation, 
and departure from thermal equilibrium However 
any produced asymmetry will be washed away by the 
SM B -|- L- violating sphaleron transitions which are ac- 
tive from temperatures 10^^ GeV down to 100 GeV [H, 
il B — L = 0. Therefore, an asymmetry vn B — L is 
generally sought which is subsequently reprocessed in 
a thermal bath via sphalerons in order to yield a net 
baryon asymmetry given by i? = a{B — L). Here, a is 
a model-dependent parameter; in case of the standard 
model (SM), a = 28/79, while in the minimal supersym- 
metric standard model (MSSM), a = 32/92 g. 

An attractive mechanism for producing B — L asym- 
metry is from the decay of the heavy right-handed (RH) 
Majorana neutrinos Since the RH neutrinos are the 
SM singlets, a Majorana mass Mjv, which violates lepton 
number is compatible with all symmetries and hence can 
be arbitrarily large beyond the electroweak scale. This 
provides an elegant way for obtaining small masses for 
the light neutrinos via the see-saw mechanism such that 
my « (mo/Mjv) [D, where TOd is the Dirac mass ob- 
tained from the Higgs vacuum expectation value (VEV). 
Moreover, a lepton asymmetry can be generated from the 
interference between the tree-level and the one-loop dia- 
grams in an out-of-cquilibrium decay of the RH neutri- 
nos, provided CP- violating phases exist in the neutrino 
Yukawa couplings. The lepton asymmetry thus obtained 
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will be partially converted to the baryon asymmetry via 
sphaleron effects. This is the standard lore for producing 
lepton asymmetry commonly known as leptogenesis . 

The present analyses of solar neutrino experiments 
favor the large mixing angle MSW solution with 

6.1 X 10-5 eV^ and tan^ 6*12 = 0.41 §, while 

3.2 X 10-3 eV^ and sin2(26l23) = (0.83 - 1) 
provides the best fit to the atmospheric neutrino data [|j . 

In addition, cosmology [|lO| , and neutrino less double- 
beta decay experiments |11[| provide an upper limit for 
the light neutrino masses. The masses and mixing an- 
gles which are required to explain solar and atmospheric 
neutrino data can be obtained in both scenarios with 
hierarchical, or quasi-degenerate neutrinos. Note that 
the hierarchical spectrum for heavy neutrinos strongly 
suggests a spectrum of light neutrinos which is hier- 
archical too unless there is a big conspiracy. On the 
other hand, a mild hierarchy of RH neutrino masses 
could be compatible with the degenerate light neutri- 
nos with a certain amount of fine-tuning. In the for- 
mer case, one may consider thermal leptogenesis sce- 
nario where heavy neutrinos come into equilibrium with 
the primordial thermal bath through Yukawa interac- 
tions. The decay of the lightest RH neutrino easily satis- 
fies the out-of-equilibrium condition by virtue of having 
a sufficiently small Yukawa coupling 0. In a model- 
independent analysis in Ref. [ p2[ , the authors have pa- 
rameterized thermal leptogenesis by four parameters; the 
CP asymmetry, the heavy RH neutrino mass, the ef- 
fective light neutrino mass, and the quadratic mean of 
the light neutrino masses. The final result was that 
an acceptable lepton asymmetry can be generated with 
Tr ~ Ah = O(lOiO) GeV, and Y.^ < V3 eV. 

However the temperature required for thermal lep- 
togenesis is marginally compatible with the maximum 
allowed one in supersymmetric theories, which is usu- 
ally constrained by thermal gravitino production |]l3|,|l^ . 
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Gravitinos with a mass ©(TeV) decay much after nucle- 
osynthesis and their decay products can change abun- 
dance of the light elements synthesized during BBN. For 
100 GeV < 7713/2 < 1 TeV, a successful nucleosynthe- 
sis requires 713/2 /s < (10^1'* - 10"n ), which translates 
into Tr < (10^ - 10^°) GcV |l|Jl|l 0. The possible ways 
for obtaining a naturally low reheat temperature include 
gravitationally suppressed decay of the inflaton in mod- 
els of high scale inflation , low scale inflationary mod- 
els [ p3|p4[p5[ | , a brief period of late thermal inflation , 
or, a completely new paradigm "reheating through the 
surface evaporation" which works even for high scale in- 
flationary models |27j . 

When the light neutrinos are almost degenerate m^^i w 
my^2 ~ Tn^^s, which requires quasi-degenerate heavy neu- 
trinos, the out-of-equilibrium condition in thermal lepto- 
genesis scenario cannot be satisfied in the minimal see- 
saw model |7|- More complicated models are required 
in this case p8|]. On the other hand, if the mass split- 
ting of the RH neutrinos becomes less than their de- 
cay widths, the perturbative calculations obviously break 
down. Then, the effect of finite decay widths of the RH 
neutrinos must be taken into account The careful 

treatment of Ref. shows that a resonant enhance- 
ment of lepton asymmetry occurs in this case, while as 
expected, it vanishes in the limit of exactly degenerate 
neutrinos. This effect can be utilized to bring down the 
scale of heavy neutrino masses, and hence the leptogen- 
esis scale po|| . 

However for almost degenerate heavy neutrinos, i.e. 
where the mass splitting is larger than the decay width, 
one has to seek non-thermal leptogenesis (which works 
for the hierarchical neutrino masses as well) in the mini- 
mal models. In this scenario RH neutrinos are produced 
non-thermally from the inflaton decay. This can occur 
during reheating if the inflaton decays to the RH neutri- 
nos, which are lighter than the inflaton, with a consid- 
erable branching ratio |31[|. Heavy neutrinos can also be 
produced via preheating |32[ (a stage of reheating where 
a resonant production of massive and/or massless bosons 
and fermions takes place p3[), or, tachyonic preheat- 
ing [Q, even if the mass of boson/fermion exceeds that 
of the inflaton. All these are rather model-dependent and 



their main features can significantly vary from model to 
model. This is the prime reason why we do not pursue 
leptogenesis via preheating mechanism here. 

In supersymmetric models one has the RH sneutrinos 
in addition. The sneutrinos are produced along with 
neutrinos during reheating, and with much higher abun- 
dances in preheating, thus serving as an additional source 
for leptogenesis Moreover, the RH sneutrinos can 

acquire a large VEV during inflation if their mass is 
less than the Hubble expansion rate during inflation Hj. 
Such a condensate starts oscillating once H{t) ~ M^r, 
therefore automatically satisfying the out-of-equilibrium 
condition. The decay of the sneutrino condensate can 
then yield the desired lepton asymmetry in the same 
fashion as neutrino decay |36| , or via Affleck-Dine mecha- 
nism |37t . This last scenario has an additional advantage 
that it solves the fine-tuning problem in _F-term hybrid 
infiationary model in a very natural way [^8|. 

The success of all these scenarios, but preheating and 
the Affleck-Dine oriented model, requires that the infla- 
ton is heavier than the RH (s)neutrinos (in the hierarchi- 
cal case inflaton only needs to be heavier than the light- 
est RH (s)neutrino.). Moreover, all the above scenarios 
are based upon the decay processes. An attractive pro- 
posal was recently made, where the lepton asymmetry 
in the visible sector is generated from the RH neutrino- 
mediated scattering of the SM Higgs and Icptons into a 
depleted hidden sector , rather than the decay of the 
on-shell heavy neutrinos. 

In this paper we propose a simple supersymmetric 
model for non-thermal leptogenesis without any need of 
preheating mechanism. In this model the inflaton is 
directly coupled to nearly degenerate RH (s)neutrinos 
which are heavier than the inflaton. Then the inflaton 
decay to the SM fields, via off-shell RH (s)neutrinos, re- 
heats the Universe and naturally leads to a sufficiently 
low reheat temperature. This same channel is also re- 
sponsible for producing the lepton asymmetry. 

In the next section we introduce our model and high- 
light several of its advantages. Then we turn to reheating 
and generation of the lepton asymmetry in this model 
and present our main results. In particular, we point 
out marked differences from leptogenesis with on-shell 
(s)neutrinos. Finally, we conclude the paper with a brief 
summary. 



* Recently, non-therrnal iimdjiction of helicity ±3/2 jl^], and 
helicity ±1/2 gravitinos [ p,6|jl7 1 from inflaton oscillations have 
been considered. For a single chiral multiplet the helicity 
±1/2 gravitino is the superpartner of the inflaton known as 
inflatino. The decay channels of inflatino have been discussed 
in Ref. J!^ | . Also, it has been suggested and explicitly 
shown |l5|f ^ that in realistic models with two chiral multiplets 
the helicity ±1/2 gravitino production is not a problem, so 
long as the inflationary scale is sufficiently higher than the 
scale of supersymmetry breaking in the hidden sector and 
the two sectors are weakly coupled. Gravitinos can also be 
produced directly from the inflaton decay |2Q|, and from the 
decay of the heavy stable neutral particles |2l| . 



II. THE MODEL 

We start by introducing our model in a supersymmet- 
ric set up. The relevant part of the superpotential is 
given by 

W D im0$$ + ig^NN ± /z,NH„L ± ^A/atNN . (1) 

Here $, N, L, and H stand for the inflaton, the RH 
neutrino, the lepton doublet, and the Higgs (which gives 
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mass to the top quark) superfields, respectively. Also, 
and M»,denote inflaton and RH (s)neutrino masses, 
respectivelytl. We assume that the inflaton is coupled 
to the RH (s)neutrinos via Yukawa coupling and h de- 
notes a typical neutrino Yukawa coupling. For simplicity, 
we have omitted all indices in h matrix and superfields, 
and work in the basis where the Majorana mass matrix is 
diagonal. Further simplifications can be made for almost 
degenerate RH (s)neutrinos where is essentially the 
same for all of them. It is also conceivable in this case 
that the inflaton is coupled with the same strength to 
three RH (s)neutrinos. This is particularly true when 
the inflaton has a nonzero VEV at the minimum which 
provides masses to the RH (s)neutrinos. We focus on su- 
perheavy RH (s)neutrinos, i.e. assuming that Mn ^ 

Now let us discuss the merits why we seek RH 
(s)neutrinos heavier than the inflaton. If Mn < m^, then 
one can easily produce on-shell (s) neutrinos from the 
inflaton decay, either perturbatively or via preheating. 
First consider (s)neutrino production in perturbative in- 
flaton decay. A perturbative decay requires a small cou- 
pling to the (s)neutrinos. This is naturally achieved when 
the inflaton lies in a hidden sector which is only gravita- 
tionally coupled to the SM sector |^^. In this case the 
total decay rate of the inflaton is given by Fd ~ m^/Mp, 
while the partial decay rate to (s)neutrinos is given by 
r^^N ~ m^Ml/M^ where Mp = 2.4 x GeV 
is the reduced Planck mass. This results in a branch- 
ing ratio ~ (Mtv/to^)^, which is too small if Mn ^ m^. 
Note that a successful leptogenesis requires an acceptable 
branching ratio given the entropy generation from reheat- 
ing, thus implying that Af^v must not be much smaller 
than m^. 

Besides, a small coupling g, which is required to en- 
sure a perturbative treatment, leads to another poten- 
tial problem. The sneutrino field N can acquire a large 
VEV during inflation. Since N is directly coupled to the 
inflaton, it might even ruin the flatness of the inflaton 
potential. On a lighter note, {N) remains non-vanishing 
after the end of inflation in any case and may contribute 
to isocurvature density perturbations | |4l[ |. This requires 
a delicate treatment of a coupled system which depends 
on the choice of a model. This is an issue which has been 
side-lined in most supersymmetric models of non-thermal 
leptogenesis, except Ref. [|^. 

If g is sufficiently large, (s)neutrinos may be produced 
in a non-perturbative manner during the stage of (ta- 
chionic) preheating p^ , p^ . For the superpotential in 
Eq. (]l|), the necessary condition for preheating reads 
gcf)!^ > m,j,, where (f)o is the initial amplitude of the in- 
flaton oscillations. This guarantees that N is heavier 



t Actually, denotes the frequency of the inflaton oscilla- 
tions around the global minimum of the potential. In models 
of chaotic inflation ~ Hi. While in new and hyhrid infla- 
tionary models it is usually (much) larger than Hi Mn. 



than the inflaton during inflation, and hence (N) = 
after the end of inflation, resulting in a simpler initial 
condition in post-inflationary era. On the other hand, 
both RH neutrinos and sneutrinos can be produced via 
preheating (sneutrinos much more abundantly by virtue 
of obeying the Bose statistics |Q). However, as men- 
tioned earlier, this is rather model-dependent. For exam- 
ple, if the inflaton has a VEV at the minimum, denoted 
as V, then it is hard to envisage an efhcient production 
of (s)neutrinos through parametric resonance. The rea- 
son is that Mn = gv and (jjQ — v in this case, which 
implies g4>o ~ M^v- It is therefore evident that there will 
be no preheating of (s)neutrinos for Mjv < m^. On the 
other hand, for Mm> preheating is possible only 
if g(l>Qm^ ^ |33|j42|] . In particular, resonant crft- 

ation rapidly ceases to be efhcient for Mn > lOm^ l42]lll. 
In the tachyonic preheating scenario, too, the produced 
(s) neutrinos usually have an abundance much less than 
the inflaton abundance when Afjv S> . In conclu- 

sion, it is very difficult (if not impossible) to obtain the 
desired lepton asymmetry in a wide range of inflationary 
models, by solely relying on non-perturbative dynamics. 

Now we count upon the advantages of our model. First 
of all, for Af^r 3> the post-inflationary dynamics is 
simpler since {N) = at the end of inflation. The Uni- 
verse is reheated through the inflaton decay to the Higgs 
and SM leptons via off-shell RH (s)neutrino. The decay 
rate, as we will see shortly, is suppressed as {rrifp/Mfq) . 
This naturally leads to an acceptably low reheat temper- 
ature when Mjq 3> m^. Furthermore, the inflaton decay 
alone is responsible for the generation of the lepton asym- 
metry. This makes the model minimal since leptogenesis 
is now directly connected with reheating. Also, the wash- 
ing out of the lepton asymmetry from thermal scattering 
of the SM leptons and Higgs is completely negligible since 
Tr « Mn. 

Our main focus will be on almost degenerate light neu- 
trinos, which can be derived naturally from almost de- 
generate RH neutrinos. An example of such a model is 
presented in Ref. |Q] , where neutrino masses and mixing 
compatible with the solar and atmospheric neutrino so- 
lutions are derived in the framework of democratic mass 
matrix. There the neutrino Yukawa matrix h is almost 
diagonal in the same basis as the Majorana mass ma- 
trix. This makes sense since when both are proportional 
to the identity matrix the light neutrinos come out to 
be exactly degenerate. Then by perturbing around this 
pattern, we can obtain a nearly degenerate texture. In 



'It has been shown in Ref. [^, that for a quadratic po- 
tential; ~ rri^cj)^, efficient resonant production of particles 
with a mass Afjv = lOm0 requires g(f>o > lO^m^. On the 
other hand, for a quartic potential ~ ^4''^, preheating of 
these particles practically disappears. Preheating in super- 
symmetric hybrid inflation model is also not efficient psl. 
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the calculations below, Mm and AM^r denote the nearly 
equal diagonal elements of the Majorana mass matrix and 
their typical differences respectively. Also h and 5h rep- 
resent the nearly equal diagonal elements of the Yukawa 
matrix and their differences respectively, while h' stands 
for the typical non-diagonal elements. It is assumed that 
AMat < Mn and h' < Sh < h. 

III. REHEATING THE UNIVERSE 

The main decay mode of the inflaton is to a four- 
body final state consisting of two Higgs/Higgsino- 
lepton/slepton particles (and their CP transforms). 
Since, we have assumed <C Mn, it is essential to 
find those diagrams which are least suppressed by pow- 
ers of Mat. These diagrams, shown in Fig. (1), which 
arise from the leading order terms in the effective super- 
potential after integrating out N, is given by 

Weff D ^m^$2 + ^g/i2$(H„L)(H„L). (2) 

We should therefore choose that part of the N propagator 
with a mass insertion, namely the part suppressed as 
1 /Mn (the other part of the propagator is proportional to 
m^/M'^). In diagrams below two opposite arrows on the 
N propagator represent this dominant part. Note that 
N propagator is proportional to (l/Af^) to the leading 
order. 

First, we evaluate the rate for inflaton decay with- 
out any specific assumption about Majorana masses and 
Yukawa couplings (except that g is diagonal and univer- 
sal, and nicf, ^ Mn)- Generically, the trajectory of the 
inflaton motion is a line on a complex 4> plane. We can 
therefore assume, without loss of generality, that only the 
real component of the inflaton has a VEV, thus treating 
the decaying inflatons as real fields. In addition, the SM 
particles are much lighter than the inflaton in the case 
under consideration (as will be confirmed by our results). 
Then the phase space factor for the four-body decay is 
readily found to be [l6 • 96 • {2TTf] ~ . The inflato n cou- 
pling to a given final state consisting of Lj (or Lj ) and 
Lfc (or Zfe), plus two (or Hu), is given by Y,i ^1m^- 
Here j and k stand for the lepton flavor. There is also 
a multiplicity factor for each final state which can be 
calculated easily. 

Given all possible weak isospin assignments, with 
flavour indices flxed, there exist a total of nine final 
states. Seven of them which consist of two fermions 
and two scalars are (1) L'jLlH^H'^, (2) l]lIh^H^, 

(3)L-LlH-Hl (4) L%hIh\, (5) LfLlHtHl (6) 

L^^L^H^H^, and (7) ^L^^H^hI. There are also two 

final states consisting of four scalars: (8) LjLf,H!^H!^ and 




b) AL = -2 



JK 




L 



c) AL = -2 



Fig. 1: Diagrams together with their CP transformed, 
for which AL = -)-2, represent the infiaton decay into 
two Higgs/Higgsino-lepton/slepton pairs at the leading 
order. 

Note that at each vertex in diagram (a), the production 

- — a - - b 

of L° (or L ) is accompanied by that of (or H^), 
and vice versa. In diagram (b), on the other hand, the 

- - a - b 

production of L° (or L ) is accompanied by that of 
(or i?Jj), and vice versa. This implies that final states in 
{1) — {4) and (7) can arise from diagram (a), while (6) 
arises only from diagram (b). On the other hand, the 
final state in (5) can arise from both diagrams. Finally, 
{8), (9) arise only from diagram (c). 

The rate for the inflaton decay to the final states in 
{1),{2) and {8) are the same and given by 
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ri = r2 



m 



^(2. [8-45,,]) X- 
^ 16 • 96 • (27r) 



Ehijhik 



(3) 



The constraint j < fc is imposed in order to avoid double- 
counting of the same final states. Note that the first num- 
ber inside parenthesis comes from the summation over all 
isospin states, while the second one represents the overall 
factor from the superposition of different contributions 
for each isospin assignment. 

Similarly, one can also evaluate the rates for the decay 
into other final states. The results are 



r3 = r4 = r9 = iri. 



While 



y (2 ■ 4) X . 

' 16 • 96 • (27r)^ 



2Mf 



(4) 



(5) 



overall rate is actually proportional to {rn^/M^). This 
is just the manifestation that these diagrams do not arise 
from the effective superpotential individually. Also, there 
exists a two-body decay channel for the inflaton, into 

HuHu{HuHu) or LL{LL), at the one- loop level. It can 
easily be derived by choosing (I/Mn) and (m^/M^) 
parts of A'^ propagators in diagram (a) and connecting 
the Hu{L), or L{Hu), lines. This channel has a much 
larger phase space factor {8tt)~ , while the dependence 
on g and h remains the same as in Fig. (1). However, the 
two-body decay rate is oc {wJ^/Mff). Thus by taking the 

one-loop factor (47r)~^ into account and for M^r > lOm^, 
it will eventually be smaller than that in Eq. (^. 

Finally, the infiaton can also decay into the SM fields 
via gravitational couplings with a decay rate Tg^-av ~ 
(z;/Mp)2(myM|), where v denotes infiaton VEV at the 
global minimum of the potential pl| . Such a decay rate 
can however be neglected compared to the four-body de- 
cay provided v <C My>. 



IV. THE LEPTON ASYMMETRY 



The total decay rate of the infiaton will be 



i=l 



(6) 



Let us come back to the case with nearly degenerate neu- 
trinos, where AAfjv < and h' < Sh < h. In this case 
each Ni (Ni) is dominantly coupled to the i-th lepton 
doublet, and the coupling is h. In consequence, off-shell 
Ni (Ni) mainly contributes to the infiaton decay to the 
final states with j — k = i. Then we can show that the 
decay rate will be given by 



(7) 



The infiaton completely decays when ~ Fd, where 
in a radiation-dominated Universe 
pO| , with being the effective number of relativis- 
tic degrees of freedom which is ~ 214 in the MSSM. 
Assuming that thermal equilibrium is achieved when 

~ Fd (which is justifiable; for the detailed discussion, 
see Refs. pSlpl), we obtain 



10-7/2 



Ml 



(8) 



Some comments are in order regarding our estimates 
of Fd and Tr . One might think that inflaton decaying 
into four scalars, the same as in diagram (b) except that 
Hu and L are replaced with and L, would occur at 
a rate only suppressed by two powers of Mjy. However, 
this is not the case since this leading order contribution 
is canceled out by that from another diagram and the 



In this section we evaluate the lepton asymmetry 
generated from the inflaton decay through diagrams in 
Fig. (1). First, we remind the readers that for the stan- 
dard case where the decay of on-shell neutrinos yields the 
lepton symmetry, one has 77L ~ "Eiei^nN./s), where 



and @ 



f{x) = y/x 



X- 1 



In 



1 + a; 



(9) 



(10) 



The first and second terms on the right-hand side of 
Eq. ( p^ correspond to the one-loop self-energy, and ver- 
tex corrections, respectively. For hierarchical N, the fol- 
lowing lower bound is found ESl 



< 



3 Mima 
8^ (Ho)' ' 



(11) 



where Mi and denote masses of the lightest heavy 
neutrino and the heaviest light neutrino respectively. 
Here {H°) = 174 sin/3 GeV is the VEV of i?„ in our vac- 
uum, with tan/3 defined as the ratio of {H°) and (H^)- 
On the other hand, a; « 1 for almost degenerate RH neu- 
trinos, and hence the self-energy contribution dominates. 
Then, it can be shown that (to the leading order) [H 



1 M 



ei ~ €2 



£3 



N 



M 



N 



47r /l2 AM, 



N {Hoy 



2m„ 



(12) 
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Now we come back to our case, where the inflaton de- 
cay via ofF-shell N (N) produces the lepton asymmetry. 
The net lepton asymmetry is generated from the inter- 
ference between diagrams in Fig (1), and the one- loop 
diagrams representing self-energy and vertex corrections 
to one of the N (N) propagators. Diagrams with one- 
loop correction to both N (N) legs are of higher order 
and will be subdominant. There are major differences 
which arise in the analysis compared to the on-shell case, 
as we note in this following discussion. To demonstrate 
these differences explicitly, we focus on self-energy and 
vertex corrections to the diagram (a) of Fig. (1), shown 
in Fig. (2). Similar arguments will go through for the 
inflaton decay through diagrams (b) and (c) in Fig. (1). 

Note that both HuL and HuL loops contribute to the 
self-energy correction, while only one of them is relevant 
in the vertex correction for a given final state. Also, 
recall that only loops with on-shell particles make a con- 
tribution to the resultant asymmetry. Thus the self- 
energy and vertex loops involving Ni actually represent 
s-channel and t-channel scattering of a Higgs-lepton or 
Higgsino-slepton pair via off-shell Ni, respectively. The 
center-of-mass energy available in these processes is at 
most equal to the inflaton mass. In consequence, the 
self-energy correction is simpiv twice as large as the ver- 
tex correction for <C Mff i. It can also be shown that 
only the mass insertion part of Ni propagator contributes 
to the generated asymmetry from self-energy correction 
of Ni . The diagram with mass insertion in Ni propagator 
will be irrelevant, exactly like the standard case p7[ . 

Important difference arises in comparison with the 
standard case is that there the center-of-mass energy 
in the two-body decay of Ni is simply determined by 
Mi. While here the energy flowing in the Ni leg is 
< -E < 7710. In the M^r limit, the Ni propagator 

is E/Mi, while the Ni propagator will simply be 1/Mi, 
see diagrams in Fig. (2). For a given final state with 
definite momenta the one-loop diagram is suppressed as 
E'^/MiMi with respect to the tree-level one. Upon per- 
forming the phase space integration over a four-body fi- 
nal state we find the suppression will be rn^/MiMi times 
some numerical factor 0(I) . For 



N 



N 



N 




a) 



HuiLk) 



'Lk{Hu) 



^T his i s similar to the a; 2> 1 limit for the standard case in 
Eq. ^ 




Lj{Hu) 



Hu{Lk) 
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- -Hu{Lk) 



b) 



Fig. 2: Diagrams representing one-loop (a) self-energy, 
and (b) vertex corrections to the decay channel shown in 
Fig. (la). The interference between these and the tree- 
level diagram results in a net lepton asymmetry. 



simplicity, we take the average energy in the Ni legs 
to be and hence the suppression comes as ~ 

nr^/AMiMi. This approximation is adequate for our pur- 
poses in the <^ Mn limit, and any difference from the 
exact result will be numerically irrelevant. The reason is 
that the main contribution to the phase space integral 
comes from the bulk of the available phase space. While, 
contribution of the parts in which the energy of some de- 
cay products is <C m^, including parts with « or 
E « m^, is suppressed. The situation will be more com- 
plicated when 7710 and Mn are not very different, since 
the energy and momentum carried by N legs is compa- 
rable to Mn- In such a case the Ni and N propagators 
can strongly depend on the phase space distribution of 
the decay products and the above approximation may 
not be sufficient. 

Now let us find the asymmetry parameter in the in- 
flaton decay. First consider diagrams in Figs, (la) and 
(2). For a given final state the tree- level and interference 
terms naturally have the same multiplicity factor. As ex- 
plained earlier the self-energy correction is twice as large 
as the vertex correction, and also, the average energy car- 
ried by each of the N legs can be approximately taken 
to be 7710/2. The contributions from both diagrams in 
Figs. (2a) and (2b) are equal, and hence the asymme- 
try receives an extra factor of 37n^/2M^, in addition to 
(I/Stt) prefactor in Eq. (||). Note that the one-loop cor- 
rection can come from each of the two N legs, which is 
equivalent to exchanging j with k. The situation will be 
similar for the asymmetry in the inflaton decay through 
diagrams (b) and (c) in Fig. (1). Thus, after summing 
over all possible final states, we obtain 



Im[(hht)„.(hht)„,(hht).;]m| 



Stt 



([hhtwi! 



(13) 



which is functionally very different from the standard 
case in Eq. (^. 

We now come back to the case with nearly degenerate 
neutrinos. Now, the denominator of the second term on 
the right-hand side of Eq. (|l|) is ~ 3h'^/Mjj. In this 
case the Yukawa matrix h is almost diagonal, and so is 
the matrix hh^. The numerator receives the main con- 
tribution from the terms with i = n and i = I td, and 
can be written as 



/i^Imffhhtl 



1 



1 



MpM, 



(14) 



This can further be sim- 

plified to 12h'^Shh'^{m'^AM n/MJ^), assuming that only 
the non-diagonal elements contain CP- violating phases. 
Since the lepton number is violated by two units in the 
inflaton decay, we finally have 



riL 3 Shh''^ AMn ( 
770 n h Mn \Mn 



(15) 



An important observation is that here the final asym- 
metry is proportional to AMn, contrary to the on-shell 
case in Eq. (|lj). Therefore the generated asymmetry 
actually decreases as the RH (s)neutrinos become more 
degenerate. This is not difficult to understand as the 
available energy in the inflaton decay tti^, is far below the 
mass of the RH (s) neutrinos Mn, independently of how 
degenerate the latter ones are. 

The total asymmetry in the baryons (after taking into 
account of sphaleron effects) can be expressed as 



nB \ 
J \ s 

1 6hh''^ AM, 
M 



N 



N 



( MnttIij 



M^ 



7770 



(16) 



where s = {2tt'^ /4:5)g^,T^. Here tu/s denotes the dilution 
from reheating. By using Eq. M), and the relationship 
rrii, ~ {h^{H^)'^ /Mn), we eventually obtain 



4.10-49/2 



5hh'^ AMn mTMp^^ 



9- 



Mn M%{H^S 



(lGeV)^ (17) 



where we have taken 777,^ « 0.1 eV. We also assume 
(i/O) = 174 GeV in below. Moreover, for AM n ^ Mn 
and as long as h' < Sh, it is sufficient to have 5h/h 



**Terms with n = I and i = n = I are real, and hence do not 
contribute to the asymmetry 
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AMn /2Mm in order to obtain degenerate light neutrino 
masses. Therefore, we may further simplify Eq. to 
find 

(18) 

Let us now present some numerical examples for nearly 
degenerate superheavy RH (s)neutrinos, i.e. Mn > 
lOm^ and AMn ^ Mn- With Mm = lOm^ and 
10~^ < h! /h < 1, the desired baryon asymmetry can 
be obtained for the range of parameters 10"'^ < <? < 1 
and 10" GeV < < 10" GeV, which result in 

10^ GeV < Tr < 10*^ GeV. With Mn = lOOm^, and 
10~^ < h'/h < 1 as before, an acceptable asymmetry 
is yielded for g = 1 and ~ 10^^ — 10^^ GeV, with 
10^ GeV < Tr < 10^ GeV. 

The merits of our model are already evident from 
these numbers. First of all, the reheat temperature is 
low (more than) enough to avoid the gravitino problem. 
Moreover, Tr <ti Mn guarantees that lepton number vio- 
lating scattering of the SM particles are completely neg- 
ligible. Especially, keeping in mind that in the MSSM, 
there are a large number of scattering processes which 
can considerably attenuate the obtained asymmetry if 
the reheat temperature Tr is close to AIn 0. In our case, 
obtaining a sufficiently low reheat temperature is more 
than welcome in this regard. Also, the robustness of the 
inflaton mass lies in a range compatible with both 
high and intermediate scale inflationary models, though 
slightly favoring high scale models, thus making the sce- 
nario more flexible. 

We shall re-emphasize the marked difference from lep- 
togenesis with on-shell (s)neutrinos, namely suppression 
of the yielded asymmetry as AMn/Mn decreases. This 
implies that our scenario works well for nearly degener- 
ate neutrinos (and perhaps even better in the hierarchi- 
cal case), while producing too little asymmetry for highly 
degenerate ones. Note that no resonant enhancement of 
the lepton asymmetry of the type discussed in Ref. ||29| ] 
will occur. However, we can expect a qualitatively simi- 
lar effect if (at least) one of the RH sneutrinos is almost 
degenerate with the inflaton. 

And the final comment before closing this section. A 
small number of on-shell (s)neutrinos might also have 
been produced non-perturbatively, from an inefficient 
preheating, and hence contribute to the resultant asym- 
metry through their decay. The asymmetry yielded in 
the decay of on-shell particles, denoted as rj^, will be 

Note that the asymmetry parameter for on-shell 
(s)neutrinos is dominated by the self-energy correction, 
given in Eq. (p^, and hence rj^ does not contain the 



suppression factor {m^/2MN) ■ On the other hand, a 
factor of 4 will be lost, relative to the off-shell case, since 
the one-particle decay of on-shell N and TV violate the 
lepton number by one unit. Thus, with AMn — Mn, 
possible contribution from on-shell (s)neutrinos can be 

neglected, provided {n^ + tin) < (Sml/M'^^n^. For 

the range of parameters considered above this is gcneri- 
cally the case. 

V. CONCLUSION 

In this paper we have provided a simple example 
for non-thermal leptogenesis with nearly degenerate su- 
perheavy RH neutrinos in a supersymmetric set up. 
We assumed that the inflaton is lighter than the RH 
(s)neutrinos, thus naturally avoiding some potential 
problems which can naturally arise. The inflaton de- 
cay via off-shell (s)neutrinos reheats the Universe and 
the model is minimal in a sense that the same channel 
is also responsible for generating the lepton asymmetry. 
As usual, the asymmetry arises from the interference be- 
tween the tree-level and the one-loop diagrams represent- 
ing self-energy and vertex corrections of (s)neutrinos, al- 
though off-shell in our case, provided neutrino Yukawa 
couplings contain CP-violating phases. However, there 
are important differences from leptogenesis with on-shell 
(s)neutrinos, which we have pronounced here. The self- 
energy and vertex corrections are now of the same order 
regardless of the degree of degeneracy. Most notably, the 
asymmetry parameter is found to be linearly proportional 
(rather than inversely in the on-shell case) to the mass 
difference of the RH (s)neutrinos. This results in a lep- 
ton asymmetry which gets smaller as the RH (s)neutrinos 
become more degenerate. 

Finally, we brieffy emphasize on remarkable advan- 
tages of this model. First of all leptogenesis can be 
accommodated rather simply without relying on non- 
perturbative production of RH (s)neutrinos. It is par- 
ticularly attractive that the desired baryon asymmetry 
can be directly generated in the final stage of reheating 
which is perturbative, regardless of any model-dependent 
effects which might have resulted in a first stage of non- 
perturbative reheating. Secondly, the suppressed decay 
of the inflaton naturally leads to an acceptably low re- 
heat temperature, which is compatible with the gravitino 
bound and also prevents any wash out of the yielded 
asymmetry. Also, with nearly degenerate (s)neutrinos, 
the desired lepton asymmetry can be generated for a 
range of inflaton mass accessible in large and interme- 
diate scale models of inflation. 

Qualitatively, we expect that this scenario also work 
(even better) in the case of hierarchical RH (s)neutrinos. 
However, a more careful study should be performed in 
order to compare the quantitative results with those ob- 
tained here. It will also be interesting to study possible 
enhancement of the lepton asymmetry when the inflaton 
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is almost degenerate with some of the RH sneutrinos. 
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